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We forecast the future constraints on scale-dependent parametrizations of galaxy bias and their 
impact on the estimate of cosmological parameters from the power spectrum of galaxies measured in 
a spectroscopic redshift survey. For the latter we assume a wide survey at relatively large redshifts, 
similar to the planned Euclid survey, as baseline for future experiments. To assess the impact of 
the bias we perform a Fisher matrix analysis and we adopt two different parametrizations of scale- 
dependent bias. The fiducial models for galaxy bias are calibrated using a mock catalogs of Ha 
emitting galaxies mimicking the expected properties of the objects that will be targeted by the 
Euclid survey. 

In our analysis we have obtained two main results. First of all, allowing for a scale-dependent 
bias does not signihcantly increase the errors on the other cosmological parameters apart from the 
rms amplitude of density fluctuations, ag, and the growth index 7 , whose uncertainties increase by 
a factor up to two, depending on the bias model adopted. Second, we find that the accuracy in 
the linear bias parameter bo can be estimated to within 1 - 2 % at various redshifts regardless of the 
hducial model. The non-linear bias parameters have significantly large errors that depend on the 
model adopted. Despite of this, in the more realistic scenarios departures from the simple linear 
bias prescription can be detected with a ~ 2 ct significance at each redshift explored. 

Finally, we use the Fisher Matrix formalism to assess the impact of assuming an incorrect bias 
model and found that the systematic errors induced on the cosmological parameters are similar or 
even larger than the statistical ones. 


I. INTRODUCTION 

In the next future large experiments like DESI |3S| and Euclid m will use galaxy clustering to obtain simultaneous 
informations on the geometry of the Universe and the growth rate of density fluctuations by measuring the galaxy 
power spectrum or the two-point correlation function at different cosmic epochs. These studies will be based on 
large surveys of extragalactic objects and will allow to accurately estimate cosmological parameters, among which 
the contribution of the Dark Energy to the cosmic density, the nature of this elusive component and, finally, to 
test non-standard theories of gravity. Since one typically observes the spatial fluctuation in the galaxy distribution, 
not in the mass, some independent phenomenological or theoretical insight of the mapping from one to the other is 
mandatory. This mapping, which is commonly referred to as galaxy bias, parametrises our ignorance on the physics 
of galaxy formation and evolution and represents perhaps the most serious source of uncertainties in the study of the 
large scale structure of the Universe. 

Galaxy bias should not regarded as a simple “nuisance” parameter. Its estimate is not just necessary to obtain 
unbiased cosmological information. However, it also allows to discriminate among competing models of galaxy forma¬ 
tion and the physical processes that regulate the evolution of stars and galaxies. Current limitations in the theoretical 
models of galaxy evolution do not allow to predict galaxy bias with an accuracy sufficient to constrain dark energy or 
modified gravity models [^. The alternative approach is phenomenological: i.e. to estimate galaxy bias directly from 
the data. 

Ref. [12], among others, has shown that future galaxy redshift surveys contain enough informations to break the 
degeneracy between the galaxy bias, clustering amplitude and the growth factor, effectively allowing to estimate galaxy 
bias from the data itself. This come at the price of making some assumption on the biasing relation. In this work we 
assume then that the bias relation between the galaxy density field and the mass field is local, deterministic but not 
necessarily scale-independent. The hypotheses of a local and deterministic bias break down on galactic scales where 
stellar feedback processes and, more in general, baryon physics become important. These scales are small compared 
to the typical size of current and future galaxy surveys and can be effectively smoothed out by focusing on large scale 
clustering analyses. 

A scale-dependent and redshift-dependent bias is a well established observational fact. Its evidence has emerged 
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from different type of clustering analyses ranging from 2 -point statistics m nn m EOl [3SM iU iZl US higher 
order statistics |18l 1241 Kill 1451 flH] . galaxy counts mugni] and gravitational lensing iHumnsiiii]. In our analysis we 
therefore assume that galaxy bias is scale and redshift dependent. In doing so we shall adopt a parametric approach in 
which one assumes a priori a theoretically justified analytic model for the scale dependent bias. In this approach we 
ignore a possible scale dependent bias deriving from primordial non-Gaussianity nmniiin] whose signature, perhaps 
also detectable in next generation surveys [19], is however supposed to be seen on scales much larger than those 
relevant for galaxy bias. 

In this work we forecast the errors on cosmological and galaxy bias parameters and assess the robustness of our pre¬ 
dictions against the choice of the bias and the fiducial model. More specifically, we consider two different parametriza- 
tions for the scale-dependent bias: a simple power-law model and the polynomial model proposed by |7j. Both provide 
a reasonable good fit to mock galaxies similar to those that will be targeted by Euclid. We also adopt two different 
fiducial bias models: a simple but rather unrealistic unbiased model, in which galaxies trace mass, and a more realistic 
model in which the bias parameters allow to match the 2-point clustering properties of the mock galaxies. Finally, we 
restrict the analysis to large scales to avoid strong non-linear effects, which allow us to consider linear prescriptions 
complemented by some additional term to account for mildly non-linearities. 

As dataset we assume a wide spectroscopic galaxy redshift survey spanning a large redshift range and consider, as 
a reference case, the upcoming Euclid survey m- 

The layout of the paper is as follows. In section |n| we provide the theoretical framework of our Fisher matrix 
analysis and discuss the bias models used. In section |III| we present the result of the analysis when galaxies are 
assumed to trace the underlying mass density field at all redshifts. In section |IV| we adopt a more realistic fiducial 
model for galaxy bias, calibrated using mock galaxy catalogs and repeat the analysis of the previous sections. The 
main results are summarised and briefly discussed in the last section. 


II. THEORETICAL SETUP 


Following |3S] we model the observed galaxy power spectrum, Pots at the generic redshift z as: 


Pobs{z,k) 


G'^{z)b'^(k,z) 

-kPshotiz) 




DAfiz)H{z) 

Dl{z)Hf{z) 


Pof{k) 


( 1 ) 


where Da is the angular-diameter distance, H{z) is the so called expansion history, i.e. the Hubble constant at 
redshift z, G{z) is the linear growth function normalized to unity at 2 = 0, f{z) = dlogG/d\oga is the growth rate, 
b{z,k) is the scale-dependent bias, Po{k) is the matter power spectrum at the present epoch, p, is the cosine angle 
between the wavenumber vector k and the line of sight direction, Pghot is the Poisson shot noise contribution to the 
power spectrum and the subscript / identifies the fiducial model. We parametrize the growth rate as / = with a 
constant growth index 7 . 

We perform the forecast using the Fisher matrix information method, i.e. by approximating the likelihood as a 
Gaussian in the parameters around a particular fiducial model, i.e. a value of the parameters that is assumed to 
approximate the 2-point clustering properties of galaxies in the real Universe. Since we are mostly interested in 
constraining galaxy bias, we rounded off the values of the cosmological parameters of the fiducial model, rather than 
using latest experimental figures. We set Hq = 0.7, ilmo = 0.25; U{,o = 0.0445; ilko = 0; primordial slope ris = 1; and 
the dark energy equation of state Wq = —0.95. Finally, we set 7 = 0.545 and rms density fluctuation at 8 h~^ Mpc 
tTg = 0.8. 

Errors in the measured spectrometric redshift, 6z « 0.001(1 -I- z), propagate into errors in the estimated distances, 
tJr = Under the hypothesis that they are independent on the local galaxy density we can model their effect on 

the estimated power spectrum as 

P (z, k) = Pobs [z, k) . (2) 

This term removes power on small scales and effectively damp nonlinear effects on wavenumber larger than . In 
order to account for mildly nonlinear effects we follow dsiiin] and multiply the power spectrum by the damping factor 


exp 



2 



( 3 ) 


This elliptical Gaussian function models the displacement field in Lagrangian space on scales ^ 10 /i ^Mpc, where 
we focus our analysis. In this expression E_l and Ey represent the displacement across and along the line of sight. 
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respectively. They are related to the growth factor G and to the growth rate / through Ej_ = EqG and Ey = 
EoG(l + /). The value of Eg is proportional to erg. For our fiducial cosmology where erg = 0.8 , we have Eg = 11 h~^ 
Mpc. This implies that power is removed below this scale. 

Assuming that the fluctuation Fourier modes are Gaussian variates, the Fisher matrix at each redshift shell is 

[II Hi 


Fij = 2 tt 



d\ogP{K)d\ogPikn) fc" 

de, d0j ' ' SttS 


(4) 


where the derivatives are evaluated at the parameter values of the fiducial model. Here, the maximum frequency 
kmax{z) is set by the scale at which fluctuations grow nonlinearly while kmin{z) by the largest scale that can be 
observed in the given redshift shell. We set a hard small-scale cut-off kmax = 0.5Mpc at all redshifts which, 
together with the damping terms ([^ and ([^ account for non linearities. On large scale we set kmin = 0.001 h~^ Mpc. 
However, its precise value is not very relevant since the contribution of low k modes to the Fisher matrix is negligible. 
Veff indicates the effective volume of the survey defined as: 


Kff = 


f 

n{r)P {k, fi) 

2 

nP (fc, fi) 

1 
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(5) 


where n = n(z) is the galaxy density at redshift z. The second equality in Equation ([^ holds if the co-moving number 
density is constant within the volume considered. This assumption, which we adopt in our analysis, is approximately 
true in a sufficiently narrow range of redshifts. For this reason, we perform the Fisher matrix analysis in different, 
non-overlapping redshift bins listed in Table |T[ together with their mean galaxy number density. The redshift range, 
the size of the bin and the number density of objects roughly match the analogous quantities that are expected in the 
Euclid spectroscopic survey. To improve the correspondency, we multiply the galaxy number densities in Table jl] by 
an "efficiency" factor 0.5 and assumes a survey area of 15,000 deg^. 

For additional robustness of our results, we marginalize over the Alcock-Paezinsky parameters. That is, we convert 
the wavenumber norm k and direction cosine /i from the fiducial to the any other cosmology using a free Hubble func¬ 
tion and angular diameter distance parameters and marginalize over them, instead of projecting over the background 
parameters flmOi NO¬ 


TABLE I: Redshift bins used in our analysis and their mean galaxy number density. Col. 1: central redshift of each redshift 
shell with width Az = 0.2. Col 2: Mean number density of objects in h^Mpc“^. These numbers match those expected for a 
Euclid-like survey according to m- 


Z '^dens 

0.6 3.56 X 10"^ 
0.8 2.42 X 10"^ 
1.0 1.81 X 10"^ 
1.2 1.44 X 10"^ 
1.4 0.99 X 10"^ 
1.6 0.55 X 10"^ 
1.8 0.29 X 10"^ 
2.0 0.15 X 10"^ 


A. Analytic models for scale-dependent bias 

The final ingredient in Eq. Q and the focus of this paper is the scale-dependent galaxy bias b{z, k). Several authors 
have proposed different models, both phenomenological and theoretical ciiiaEsinTiiMiiig. In this work we are not 
too concerned on the accuracy of bias models. Our goal is to assess the impact of a scale-dependent galaxy bias in 
the analysis of future galaxy surveys. For this purpose we have decided to adopt two rather simple models, the Power 
Law and the Q-Model, that nonetheless provide a good match to the galaxy bias measured in numerical experiments, 
as we shall see. The reason for choosing these models is twofold. First, they have been already used in the literature, 
making it possible to compare our results with existing ones and use previous results to set the range in which the 
model parameters can vary. Second, their simple form allows us to compute the power spectrum derivatives in the 
Fisher matrix analytically. 
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TABLE II: Bias parameters for Type 1 fiducial models. 




FMl-PL 

FMl-Q 


n = 0 

n = 1 n = 2 



bo 

bo bi bo bi 

bo Q A 

all 

1 

10 10 

1 0 0 


The Power Law bias model has the form HZ]: 

b{z,k) = bo{z)+ bi{z) \^] , 


( 6 ) 


where the pivot scale ki is introduced only to deal with dimensionless parameters. Its value does not impact on our 
analysis and, without lack of generality, we set ki = 1 h Mpc“^. The slope n is not treated as a free parameter but is 
kept fixed. However, to check the sensitivity of our results on the power law index we have considered three different 
values: n = 1, 1.28, 2. As we shall see the value n = 1.28 corresponds to the one that provides the best fit to the bias 
of mock galaxies measured in simulated catalogs, n = 1 also provides an acceptable fit to the mock galaxy bias. The 
case n = 2 should be regarded as an extreme case since it provides a poor fit to the simulated data both at large and 
at small scales. We note that the Power Law model is similar to the one proposed by jSS] in those fc—ranges in which 
the power spectrum can be approximated by a power law. 

The Q-Model is also phenomenological. It has been proposed by (7) from the analysis of mock halo and galaxy 
catalogs extracted from the Hubble volume simulation. Its form is 


b{z,k) = bo{z) 


1 + Q{z){k/kiy 


1 + A{z){k/ki) _ 


1/2 


(7) 


In our analysis all three parameters 6o, Q and A are free to vary in each redshift bin. Therefore, the Q-Model has 
additional degrees of freedom with respect to the Power Law model. 

Finally, we need to specify the parameters of the fiducial model. In this work we consider two different fiducial 
models. Type 1 and Type 2, corresponding to two different choices of parameters for each bias model, totalling to 
four fiducial models. Type 1 models (denoted as FMl-PL for the power law bias and FMl-Q for the Q-model) 
represent the simple but rather unphysical case of galaxy tracing mass at all redshifts. Let us stress that assuming a 
scale-independent fiducial model does not imply that we are assuming a scale-independent bias since we differentiate 
the power spectrum in the Fisher matrix with respect to the scale-dependent bias coefficients: bi, Q, A. Instead, 
assuming a scale-independent fiducial model only implies that in Type 1 models the derivative is evaluated at the 
fiducial values &i, QandA — 0. The parameters that identify the fiducial models are listed in Table |ll| Note that we 
also consider for comparison the case of scale-independent bias (first row). This case is identical to choosing n = 0 in 
the power law model, i.e. to b{z) = bQ{z) -|- bi{z), so we will refer to this case as n = 0 fiducial. 

Type 2 models (denoted as FM2-PL and FM2-Q) are more realistic. The fiducial parameters of the Power Law 
and Q-Models were determined by matching the bias of mock Ha emitting galaxies in the simulations described in 
Section jlV A| Table jlHj lists the parameters of the fiducial models. 

Note that for the power law cases the values of the parameters depend on the choice of the power law index n so 
that, since we explore the cases n = 1,1.28, 2, we end up by having several Power Law fiducial models. 


TABLE III: Bias parameters for Type 2 fiducial models. 





FM2-PL 



FM2- 

Q 


n = 

= 1 

n = 

1.28 

n - 

= 2 




2 : 

bo 

bi 

bo 

bi 

bo 

6 i 

bo 

A 

Q 

0.8 

1.04 

0.67 

1.09 

0.66 

1.17 

0.68 

1.26 

1.7 

4.54 

1.0 

1.13 

0.74 

1.19 

0.75 

1.28 

0.79 

1.36 

1.7 

4.92 

1.2 

1.22 

0.99 

1.30 

0.97 

1.41 

1.02 

1.49 

1.7 

5.50 

1.4 

1.36 

1.09 

1.44 

1.06 

1.55 

1.12 

1.63 

1.7 

5.70 

1.6 

1.49 

1.22 

1.58 

1.19 

1.71 

1.25 

1.75 

1.7 

6.62 

1.8 

1.61 

1.40 

1.72 

1.40 

1.88 

1.44 

1.92 

1.7 

6.99 
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The complete set of parameters that are free to vary is then /i, OhO: ^mOj 7 j plus, for each redshift shell, 

the shot noise Pnoise and bo, bi or alternatively bo, A, Q. This amounts to a total of 39 independent parameters when 
we adopt 8 redshift bins. All the other parameters (the bias power law index n, the damping model, the survey 
parameters like galaxy density, volume, redshift error as well as the dark energy parameters wq and wi are kept 
fixed). 


B. Derivatives in the Fisher Matrix analysis 


In order to evaluate the Fisher matrix we compute the derivatives of the power spectrum in Eq. ( [T|) with respect 
to the parameters on the fiducial model. The Itr error for each parameter of the model, pi, is ap. = where 

F~^ is the inverse Fisher matrix. 

Focusing on the free parameters of the biasing function we have for the Power Law bias model: 


d\nP 

dbo 

d\nP 

dbi 


1 / 

1 / 


For the Q-model we have instead : 


2 2fp^ 

bl{bl + fp^) 
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dlnP 


dlnP 
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l + Qf{k/kiY 
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1 + Af{k/ki) 


fp^{k/kif 


l + Qf{k/k,r 


l+Q/(fc/fci)^ 

l+Af(k/ki) 
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fp^{k/kx)bl 


[l + Af{k/ki)\ 
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i+Qf(k/kiy 

l+Af{k/ki) 
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( 8 ) 

(9) 


( 10 ) 


( 11 ) 


( 12 ) 


III. TYPE 1 FIDUCIAL MODELS. RESULTS 

In this Section we show the results of the Fisher Matrix analysis performed using Type 1 models FMl-PL and 
FMl-Q. They represent the case of a survey of objects that trace the underlying mass density field in an unbiased 
way. 


A. Power law case 

For the power law model we have explored two cases corresponding to different choices of the power law index: 
n = 1 and n = 2. The more realistic case n = 1.28 will be considered in the next Section. In addition, we consider the 
case n = 0 that represents the scale-independent hypothesis often assumed in many Fisher matrix forecast papers. 
For this particular case we do not compute derivatives with respect to bi. The parameters of the fiducial models are 
reported in Table 

The goal of our analysis is twofold: to assess the impact of a scale-dependent bias on the measurement of the 
cosmological parameters and to estimate the accuracy with which we can measure the parameters that characterise 
the bias. Let us first focus on the first task. 

The Icr errors on the cosmological parameters are listed in Table |IV| Unless otherwise specified, the quoted errors 
are always obtained after marginalising over all the other parameters. For all parameters except the mass variance erg 
and the growth index 7 the errors are largely independent from n. In fact in most cases they slightly decrease when 
the scale dependence is stronger. The values of erg and 7 show the opposite trend, although the effect is quite small 
(below 10 %). We conclude that allowing for a scale dependent bias has little effect on the precision in which we can 
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FIG. 1: 68 % probability contours for erg and 7 . Black, continuous: standard scale-independent case, i.e. n = 0. Red, Dotted: 
FMl-PL with n = 1. Cyan, continuous: FMl-Q. Blue, Dot-Dashed: FM2-PL with n — 1. Green, dashed: FM2-Q. 


measure most cosmological parameters. It is interesting to note that the accuracy of the growth rate 7 is 4-5% when 
marginalising over all parameters, including the scale and redshift-dependent bias. Figure [l] gives a visual impression 
of this fact. It shows the 1-a likelihood contours for erg and 7 for most of the fiducial models explored in this paper. 
The likelihood ellipse obtained in the FMl-PL case when n = 1 (red, dotted curve) is only slightly larger than that 
corresponding to a scale-independent bias (continuous, black line). The likelihood contours obtained for n = 2, not 
plotted to avoid overcrowding, are also similar. We further notice that there is little correlation between erg and 7 . 


TABLE IV: la errors on cosmological parameters for Type 1 fiducial models. 


Error 

n = 0 

FMl-PL 

n = 1 n = 2 

FMl-Q 

ah 

0.036 

0.038 0.037 

0.039 


0.015 

0.016 0.015 

0.016 


0.0034 

0.0036 0.0034 

0.0036 


0.036 

0.042 0.036 

0.044 

(7 'Y 

0.024 

0.025 0.028 

0.029 

<^(T8 

0.0036 

0.0044 0.0045 

0.0047 


The expected Icr errors for the bias parameters bo and bi are listed in Table [Vl for the cases n = 0, 1 and 2. Errors 
on bi are larger than those on bg and their size increase with n. When n = 2 they are ten times larger than with 
n = 1. On the contrary, errors on bg weakly depend on n. This is not surprising since bi is constrained by the power 
spectrum behaviour at high fc, the larger the value of n the larger the values of k, where our analysis is less sensitive 
due to the damping terms and the hard kmax cut. A second trend is with the redshifts: errors on the bias parameters 
increase with the redshift, irrespective of the n value. Again, this is not surprising since it merely reflects the fact that 
the effective volume of the survey monotonically decreases when moving to high redshifts due to the smaller galaxy 
densities. 

When compared to the results of Ref. [TT], in which bias was assumed to be scale-independent, we notice that our 
constraints on bo are twice weaker than their "optimistic, internal bias" case at z = 1.8 and z = 2.0. This quantifies 
the effect of allowing for an additional degree of freedom, the scale dependent bias, represented by the new parameter 
hi. 

In Figure we show the 68 % probability contours in the hg-hi plane for n = 1 and n = 2, respectively. Larger 
ellipses refer to higher redshift bins. For the case n = 1 there is a strong anti-correlation between hg and hi which 
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TABLE V: Errors on bias parameters for Type 1 fiducial models. 





FMl-PL 


FMl-Q 


n = 0 

n = 

1 

n = 

2 




2 

O'bo 

Ubg 


O'bQ 


O'bo 


CJ-A 

0.6 

0.007 

0.013 

0.14 

0.0081 

1.2 

0.017 

3.04 

0.35 

0.8 

0.008 

0.013 

0.13 

0.0093 

0.97 

0.017 

2.5 

0.32 

1.0 

0.009 

0.013 

0.12 

0.011 

0.86 

0.017 

2.2 

0.31 

1.2 

0.010 

0.014 

0.12 

0.012 

0.82 

0.018 

2.2 

0.31 

1.4 

0.011 

0.014 

0.13 

0.013 

0.91 

0.019 

2.5 

0.34 

1.6 

0.012 

0.016 

0.16 

0.014 

1.2 

0.023 

3.4 

0.42 

1.8 

0.014 

0.019 

0.22 

0.016 

1.9 

0.027 

5.4 

0.59 

2.0 

0.018 

0.026 

0.34 

0.019 

3.3 

0.037 

9.4 

0.97 


stems from the fact that an increase in the linear bias term bo can be partially compensated by reducing the amplitude 
of the scale-dependent term bi. Increasing the scale dependency, i.e. setting n = 2 reduces the correlation between 
bo and bi. This is due to the fact that a strong scale dependent bias has little impact on large (fc ^ fci) scales and 
therefore cannot effectively compensate a variation of the linear bias on the scales that are relevant for our analysis. 



FIG. 2: 68% probability contours for the parameters bo and for of the FMl-PL model. The dotted red line is the case with 

n = 1 while the dashed line in blue shows the case n = 2. The redshift bins are 2 = 0.6,1.8, 2.0, from inside out. 


B. Q-model 

This version of Type 1 fiducial model is characterised by three parameters, bo, A and Q, rather than two. We have 
repeated the same analysis performed as in the FMl-PL case and summarised the results in Tables |IV| and [V| 

The errors on the cosmological parameters are very similar to those obtained in the FMl-PL case, confirming that 
the accuracy in the estimate of the cosmological parameters is little affected by the adoption of a scale-dependent 
bias model, even when we introduce an additional degree of freedom. For this reason the 68% error contours for (jg 
and 7 (cyan thick curve in Figureare quite close to that of the FMl-PL n = 1 case (red, dotted curve). 

Errors on the linear bias parameters bo (Table are small, with a magnitude similar to that of the FMl-PL case 
with n = 1. On the contrary, the errors on A and Q are quite large, although we cannot directly compare their size 
to the errors on bi. This is not entirely surprising: it is the effect of having one more parameter to marginalize over. 
To further investigate the possible degeneracy among the bias parameters we plot the 68 % uncertainty contours in 














the A-Q plane in Fig. The size of the errors, and consequently the area of the corresponding ellipse increases with 
the redshift. They are positively correlated and the strength of the correlation also increases with the redshift. 



FIG. 3: Contours plots for 68 % probability contours for the parameters A and Q of the FMl-Q model. We decided to plot 
only the bins 2 = 0.6, 0.8,1.6,1.8, 2.0, from inside out, to improve clarity. 


IV. TYPE 2 FIDUCIAL MODELS. RESULTS 


In this Section we repeat the Fisher matrix analysis performed in the previous Section using the same models for 
the power spectrum but assuming a Type 2 fiducial model for the bias. The parameters of the fiducial models have 
been provided by the best fit to the bias of the simulated Ha galaxies that will be targeted by the Euclid spectroscopic 
redshift survey. According to Ref. m Euclid will produce a redshift catalog of Ha-line emitting galaxies with line 
flux above Hamin = 3-10“^® erg cm“^ s“^ with high completeness and purity. Ideally, one would fix the parameters of 
the fiducial model by matching observations. Unfortunately, the samples currently available are too small to reliably 
estimate the bias of high redshift, Ha galaxies. For this reason we consider instead a catalog of simulated Ha galaxies 
designed to mimic the predicted properties of the objects that will be observed by Euclid and measure their bias from 
the clustering properties. The detailed procedure is described in the next Section. 

Another important difference with Section IV is that we now consider a smaller redshift range {z = [0.8,1.8]). This 
is a conservative choice mainly driven by the goal of matching as close as possible the range that will be best exploited 
for scientific analyses in the Euclid survey. 


A. Galaxy bias from mock catalogs 

We consider the so-called "100 deg^ light-cone" mock galaxy catalog [50] obtained by applying the semi-analytic 
galaxy formation model GALFORM U to the outputs of the Millennium V-body simulation jUj. The simulated volume 
consists of a light-cone with a sky coverage of 100.206 deg^ spanning the redshift range [0.0, 2.0]. Mock galaxies 
are characterised by angular position, redshift, co-moving distances and a number of observational properties among 
which the luminosity of the Ha line. From these mocks we selected a subsample of objects in the range z = [0.6,2.0] 
and with Ha line flux larger than Hamin- 

To estimate galaxy bias we have compared the power spectrum of the Ha galaxies to that of the linear mass power 
spectrum obtained from GAME [^ in six redshift bin of width Az = 0.2. The galaxy bias and its dependence on k is 
obtained from the ratio of these two quantities. This bias estimate is quite noisy due to the limited number density 
of mock galaxies. To regularise its behaviour we fit the Power Law and Q-model and determine the parameters of the 
fiducial by minimising the difference between model and estimated bias. Dealing with analytic expression for the 
bias also facilitates its implementation in the Fisher matrix analysis. 
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The detailed procedure to determine the realistic model the galaxy bias is as follows: 


• We extract six, partially overlapping cubic boxes from the simulated light cone. These boxes are fully contained 
in the light cone, aligned along the line of sight and centred at the redshifts indicated in Table |III| Their size 
increases with the redshift. 


• We select galaxies brighter than Hamin in each box and compute a statistical weight to account for the selection 
implied by the flux cut. 

• We use a Fast Fourier Transform based estimator, similar to mi, to measure the power spectrum of the 
mock galaxies within each cube. We ignore the effect of peculiar velocities and limit our analysis between 
fcmax = l/iMpc~^ and fcmin- The value of fc^ax is chosen to minimise the impact of shot noise and aliasing. The 
value of fcniin depends on the redshift and is determined by the size of each box. 

• We estimate the scale dependent galaxy bias from the ratio between the mock galaxy power spectrum and 
the CAMB power spectrum of the mass obtained using the same cosmological parameters as the Millennium 
simulation. 


• We fit the power law model and the Q-model to the estimated mock galaxy bias, b{k, z), by minimising the 
of the residuals. In the procedure we assume that the errors on the estimated power spectrum are the same as 
of the power spectrum estimator computed in m- In the minimisation procedure we fixed some of the bias 
parameters. For the Power Law model we set n = 1, 1.28 and 2. For the Q-model we set A = 1.7 according to Ref. 
[7]. The minimisation is first carried out over a conservative range of wave numbers, [fcmin, ^max] = [0.03,0.3], 
that is gradually increased in both direction and stop when y^/d.o.f. ^ 1. Note that the value n = 1.28 is the 
one that provides the best fit to the measured bias when n is also free to vary. 


Ill 


The parameters of the fiducial models obtained from the fit listed in Tab. 

Law model the best fit to the data is obtained with n = 1.28. 
whereas n = 2 provides a poor fit both at small and large fc-values. For this reason the case n = 2 should be regarded 
as an extreme and somewhat unrealistic case. 


We reiterate that for the Power 
The quality of the fit is still acceptable for n = 1 


B. Power Law case 

Table |VI| shows the la errors on the cosmological parameters when one assumes FM2-PL as a fiducial model for 
galaxy bias. Errors are similar to those of the Type 1 model for all parameters except for erg and 7 They show 
little or no dependence on the power law index n. Errors on erg and 7 are significantly larger since these parameters 
are partially degenerate with the bias. Indeed, the amplitude of the power spectrum is sensitive to the combination 
a%b{k,z) and the strength of the redshift distortion to f/b{k,z). In addition erg and 7 are significantly anti-correlated 
(solid ellipse in Fig.[^. This anti-correlation reflects the fact that an increase in amplitude of the mass power spectrum 
(i.e. an increase in ag) can be compensated by a decrease in the growth rate The reason why this anti-correlation 
was not seen in the Type 1 model case is that in this new fiducial model the bias is significantly larger than unity. As 
a consequence a much larger increase in the growth rate is now required to compensate a variation in erg. 


TABLE VI: la errors on cosmological parameters for Type 2 fiducial models. 


Error 

FM2-PL 

n = 1 n = 1.28 n = 2 

FM2-Q 

ah 

0.035 0.035 0.029 

0.031 


0.015 0.014 0.012 

0.013 


0.0033 0.032 0.0028 

0.0027 


0.0402 0.038 0.029 

0.034 

(7 'Y 

0.036 0.035 0.036 

0.052 

CJ'fTg 

0.0049 0.0048 0.0046 

0.0073 


The errors on the bias parameters bo and bi are listed in Table [VTI| for the three power law indices considered. The 
trend with n is the same as in the Type 1 case. Errors on ho do not significantly depend on n. They do increase with 
the redshift (although the relative error abgjbo does not). Errors on bi significantly increase with n. 
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TABLE VII: Errors on bias parameters for Type 2 fiducial models. 





FM2-PL 



FM2-Q 


n = 

1 

n = 1.28 

n = 

2 




z 

C^bQ 

^bi 

^bi 

abg 

^ bi 

O'bQ 


O-A 

0.8 

0.014 

0.12 

0.012 0.19 

0.011 

0.91 

0.029 

2.6 

0.46 

1.0 

0.016 

0.12 

0.014 0.18 

0.013 

0.76 

0.033 

1.9 

0.39 

1.2 

0.018 

0.12 

0.016 0.18 

0.016 

0.66 

0.039 

1.6 

0.36 

1.4 

0.021 

0.13 

0.019 0.19 

0.019 

0.69 

0.045 

1.6 

0.34 

1.6 

0.024 

0.16 

0.022 0.24 

0.023 

0.89 

0.049 

1.8 

0.37 

1.8 

0.028 

0.21 

0.026 0.32 

0.027 

1.3 

0.056 

2.1 

0.40 


Fig.a illustrates the covariance between bg and bi. The different ellipses are centred at the best fit values bg and 
bi listed in Table |III| As their magnitude steadily increases with the redshift the center of the corresponding ellipses 
shifts accordingly. 

Errors on bg are small and show a weak dependence on the redshift and on n. Errors on hi are comparatively 
larger, show a weak dependence on the redshift and a strong dependence on n. One consequence of this is that the 
possibility to detect a scale dependent bias depend on the bias model itself. For n = 1 the scale dependency, i.e. the 
fact that hi ^ 0, is detected with a significance of more than 4cr in each redshift bin and to more that 3cr for n = 1.28, 
i.e. with the most realistic fiducial model. Instead, for the extreme case n = 2 the size of the errors hides the scale 
dependency of the bias. 


C. Q-model 

For the FM2-Q case the errors on the cosmological parameters, listed in Table |VI[ are remarkably similar to those 
of the FM2-PL case except for erg and 7 whose uncertainties are significantly larger. This is clearly illustrated by the 
68 % error ellipse in Fig. (green, dashed curve) that also confirms the anti-correlation between the two parameters 
noticed in the Type 1 analysis. In addition, as in the Type 1 case, the increase in the error size reflects the fact that 
the number of free parameters in the Q-model is larger than in the Power Law model. 

The errors on 60 , A and Q are listed in the last three columns of Table |VII[ Relative errors on bg are similar to 
those of the FMI-Q case and larger that those of the FM2-PL case, again reflecting the increase in the number of 
free parameter in the bias model. Absolute errors on A and Q are smaller than in the FMl-Q case and appear to 
be correlated, as shown in Figure]^ which is expected since the bias in the Q-model depends on their ratio ([^. The 
errors on A, Q weakly depend on the redshift whereas the value of Q significantly increases both at small and at large 
redshift (the value of A is set equal to 1.7). What is most remarkable is the fact that in this realistic fiducial model 
the bias parameter Q is significantly (^ 3 (t) different from zero, meaning that the scale dependence in galaxy bias 
can be detected even adopting a 3-parameter model. For the parameter A this significance is smaller and depends on 
the redshift: it increases from ~ Icr at z = 0.8 to ~ 2a at z = 1 . 8 . 


V. SYSTEMATIC BIAS 

Choosing a particular bias fiducial when the galaxy sample is characterized by another bias model leads to a change 
in the best fit of the other parameters. It is interesting then to estimate this “systematic bias” due to the galaxy bias 
in one of our cases. For simplicity, we perform this test only for our reference case, the power-law Type I bias, i.e. 
FMl-PL n = 1. 

Suppose then we do not estimate the bias from the galaxy clustering using our method and simply assume the fiducial 
model FMl-PL, i.e. we fix the bias at every redshift. This of course produces narrower errors on the cosmological 
parameters, but the price to pay is that the results will depend sensitively on the bias choice. Indeed, if the chosen 
bias model is wrong, all the cosmological parameters will shift by a certain amount. We wish to estimate this shift. 

The maximum likelihood estimator 9c,, given a likelihood function L[6a) for the parameters 9a, is obtained by 
solving the system of equations 


Aa =0 


(13) 
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FIG. 4: 68 % probability contours for the parameters bo and bi of the FM2-PL model for different values of n. They are centred 
at the bo and bi values indicated in Tab. m Ellipses of increasing redshifts are centred at increasing values of bo- In the left 
upper panel we show the case n = 1. In the right panel we show the n = 1.28 case. In the low left panel we show the n = 2 
case. Only the redshifts 2 = 0.8,1.2,1.6,1.8 are shown to avoid overcrowding. 


where £ = — logL. In the Fisher matrix approximation one has £ = ^DiDjPij — where P = C~^ is the precision 
matrix (the inverse of the data correlation matrix) and = di — £ is the vector of data minus theory points. 

Suppose now L depends also on a parameter s (that we refer to as a systematic parameter) that has been assigned 
a particular value, and we want to estimate how 0^ changes when s is shifted to another value s + Ss, where Ss is 
assumed very small. Then we have 


£{s + Ss) = £{s) + £.sSs 


(14) 


and the equations that give the new maximum likelihood estimator vector 9 — 9 + S9 become to first order in Ss and 
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FIG. 5: 68 % probability contours for the parameters A and Q of the FM2-Q model. We plot only the redshifts 2 ; = 
0.8,1.2,1.6,1.8, left to right, to avoid overlaps. 


S0 


C{s + Ss)^a = F,a \g -\-{C^as |e)<5s 

(15) 

= F^a \g F{C^ap \e)^0p + {F^as le)<5s 

(16) 

= [F^ap \e)^0p + {F^as le)*^® = 0 

(17) 

where we employed Eq. (13). Finally we obtain 

59a = -FjpC^psSs 

(18) 


(sum over /3). If there are several systematic parameters, then a sum over the s parameters is understood. Now, once 
we average over many data realization, {C^a/s) = F^is is the Fisher matrix, and (/I,/ 3 s) = Fps is a sort of systematic 
Fisher matrix. In practice, this means that one includes the systematic parameters s/ in a general Fisher matrix that 
contains also the parameters 0a, and then selects the a/3 and the (3si submatrices and produces the sum over /3 and 
Si, namely 11251 1 

S0C. = -FjF/ss.Ss, (19) 


We consider separately the case in which one shifts the bo and bi parameters. We consider separately a systematic 
change in bo or bi of 6si =0.1 in every bin. The results obviously scale with 6si. The larger the coefficients F~^Fps. 
are, the more sensitive a parameter is to the determination of the fiducial bias. We see from Table |VIII| that the 
largest effect occurs on Ug due to a wrong bo, a shift much larger than the statistical error (compare with Table IVI. 
In general, a wrong choice of bo of magnitude 0.1 yields a shift on the cosmological parameters from 2 to 5 times their 
statistical errors. This is a powerful argument to insert the bias parameters in the Fisher matrix (and in the data 
likelihood when real data will be available) rather than freezing them to some fiducial value. 

We notice that the size of systematic errors is larger for those cosmological parameters that quantify the amplitude 
or the overall shape of the power spectrum and comparatively smaller for those related to sharp spectral features 
like ilb for BAOs. We conjecture this will apply also to the characteristic step at the free-streaming scale induced by 
massive neutrinos. Finally, a wrong choice of bi (by the same amount of 0.1 in every bin), on the other hand, induces 
a shift quite smaller than the statistical error. Allowing for a scale dependent bias makes in this case little difference. 


VI. CONCLUSIONS 

In this paper we have investigated the impact of a scale-dependent galaxy bias on the results of the clustering 
analysis performed in next generation surveys. We focused on the galaxy power spectrum and, as a case study, we 
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bo 

61 

h 

- 0.12 

0.0016 

9.mh^ 

-0.052 

0.00088 

Q,hh? 

- 0.011 

0.000098 

ris 

0.14 

-0.013 

7 

- 0.12 

- 0.011 

Cs 

0.025 

-0.0046 


TABLE VIII: Shift of the maximum likelihood parameter estimates due to a systematic error of 0.1 in the bias parameters bo 
and 6 i in every bin. 


have considered a spectroscopic redshift survey with characteristics similar to those expected for the Euclid survey. 
In particular we focused on two issues: ij the impact of a scale dependent bias on the estimate of the cosmological 
parameters and iij the accuracy with which it will be possible to detect such scale dependence. 

For this purpose we have considered two different models for the scale-dependent galaxy bias and obtained a Fisher 
matrix forecast assuming two sets of fiducial models. Both scale dependent bias models are characterised by a linear 
bias parameter bo and additional parameters that quantify scale dependence. However, in the Power Law model one 
of them, the power law index, is fixed, whereas in the Q model both parameters, A and Q, are free too vary. All bias 
parameters change with redshift, so they should really be regarded as functions specified at different redshifts. For 
as concerns the fiducial models, we have considered two scenarios. The first one is the simple but unrealistic case of 
a population of unbiased mass tracers. The second one is that of a more realistic galaxy bias relation calibrated on 
simulated catalogs of Ha-line emitting objects mimicking the expected characteristics of the Euclid survey. 

Our main results can be summarised as follows: 

• Allowing for a scale dependent bias does not increase significantly the errors on cosmological parameters, except 
for the growth index 7 and the rms density fluctuation erg- 

More specifically, in our analysis we find that errors on ho, fimOt ^fcOj Ug are insensitive to a scale 
dependency in galaxy bias, to the specific form of scale-dependent bias and to the choice of the fiducial model. 
In fact these errors are only slightly larger than those expected when one assumes that bias is scale independent. 
On the contrary errors on 7 and erg are rather sensitive to the bias model and to the fiducial. In the ideal scenario 
of a population of unbiased mass tracers the expected relative error on 7 is smaller than 5 %. However, when a 
more realistic bias model and survey setup are adopted, the relative error increases to ^ 6 % for a 2 -parameter 
bias model and to ^ 9 % for a 3-parameter model. Clearly, an accurate estimate of the growth index would 
benefit from a reliable model for galaxy bias that can be captured by few free parameters. 

In addition, 7 and erg are correlated. This correlation is expected since the clustering analysis constrains the 
amplitude of the power spectrum, proportional to the product b{k,z)a^{z) and its redshift distortions, which 
are proportional to f{z)/b{k, z). The degree of degeneracy, however, depends on the fiducial, being larger when 
more realistic bias models are adopted. 

• The linear bias parameter bo can be determined within a few %. The relative error is rather insensitive to the 
choice of the fiducial and slightly increases with the redshift. As expected, the errors increase with the number 
of free parameters in the model and therefore is larger in the Q Model than in the Power Law one. Yet, even 
in the worst case (the realistic case of a Type 2 Q-Model bias at z = 1.8) the relative error is of the order of 
2 %. This is of the same order as the uncertainty in the estimate of bo obtained when one assumes no scale 
dependency, i.e. when one does not consider derivatives with respect to bias parameters other than bo in the 
Fisher matrix analysis (see e.g.|12|L 

• The accuracy with which one can estimate the bias parameters that describe the scale dependency depends on 
the bias model and on the fiducial. The case in which the scale dependence is stronger, i.e. the Power Law 
model with n = 2, provides the worst results. This is not surprising since in this case the scale dependence is 
pushed at small scales where our Fisher matrix analysis, optimised to probe linear to mildly nonlinear scales, 
is less sensitive. However, we should stress that this n = 2 case is hardly realistic. In fact this model provides 
a poor fit to the bias of Ha galaxies in the simulated Euclid catalog. On the contrary, the Q-model, the Power 
Law model with n = 1.28 and, to a lesser extent, the Power Law model with n = 1 match the bias estimated 
from the mocks. 

When we restrict our analysis to the more realistic cases we find that scale dependency can be clearly detected 
at all redshifts. The significance of the detection, quantified by the departures of the parameters &i, A and Q 
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from zero, depends on both the bias and the fiducial. For the Power Law model with n = 1 is larger than 4cr 
at each redshift, decreasing to ^ Scr for n = 1.28. For the Q-model scale dependence is best detected through 
the parameter Q Scr) whereas the significance of a non zero A-value is expected to be smaller and redshift 
dependent. 

• An additional benefit in including galaxy bias in the likelihood is that of reducing the impact of systematic 
errors. Using a simple argument based on the “systematic Fisher matrix” approach (arXiv:1101.1521) we have 
shown that adopting an incorrect bias model induces significant systematic errors that, for some cosmological 
parameter, can exceed statistical uncertainties. 

All in all we conclude that a significant detection of a scale dependent bias is within reach of next generation redshift 
surveys and that contemplating such scale dependence does not significantly decrease the accuracy in the estimate 
of most cosmological parameters. We note that a scale dependent bias can also be induced by primordial non- 
Gaussianity. Its signature is expected on much larger scales than those affected by galaxy bias and there should be 
ample possibility to disentangle the two effects. Yet, the need to allow for additional free parameters to account for 
non-Gaussian features may have an impact on the precision of the analysis. We shall explore this issue in a future 
paper. 


VII. ACKNOWLEDGEMENTS 

The research leading to these results has received funding from the European Research Gouncil under the Eu¬ 
ropean Gommunity’ s Seventh Framework Programme {FP7/2007 — 2013 Grant Agreement no. 279954), therefore 
EM acknowledges the financial support provided by the European Research Gouncil. EB acknowledges the finan¬ 
cial support provided by NFN-PD51 INDARK, MIUR PRIN 2011 ’The dark Universe and the cosmic evolution of 
baryons: from current surveys to Euclid’ and Agenzia Spaziale Italiana for financial support from the agreement 
ASI/INAF/I/023/12/0. LA acknowledges financial support provided by DFG TRR33 “The Dark Universe”. LA ac¬ 
knowledges useful discussions within the Euclid Theory Working Group. We wish to thank Alex Merson for providing 
the 100 deg^ Euclid light cone mock catalogs. 


[1] L. Amendola, S. Appleby, D. Bacon, T. Baker, M. Baldi, N. Bartolo, A. Blanchard, C. Bonvin, S. Borgani, E. Branchini, 

C. Burrage, S. Camera, C. Carbone, L. Casarini, M. Cropper, C. de Rham, C. Di Porto, A. Ealet, P. G. Ferreira, F. Finelli, 
,J. Garcia-Bellido, T. Giannantonio, L. Guzzo, A. Heavens, L. Heisenberg, C. Heymans, H. Hoekstra, L. Hollenstein, 
R. Holmes, O. Horst, K. Jahnke, T. D. Kitching, T. Koivisto, M. Kunz, G. La Vacca, M. March, E. Majerotto, K. Markovic, 

D. Marsh, F. Marulli, R. Massey, Y. Mellier, D. F. Mota, N. Nunes, W. Percival, V. Pettorino, G. Porciani, G. Quercellini, 
J. Read, M. Rinaldi, D. Sapone, R. Scaramella, C. Skordis, F. Simpson, A. Taylor, S. Thomas, R. Trotta, L. Verde, 

F. Vernizzi, A. Vollmer, Y. Wang, J. Weller, and T. Zlosnik. Gosmology and Fundamental Physics with the Euclid 
Satellite. Living Reviews in Relativity, 16:6, September 2013. 

[2] P. Arnalte-Mur, V. J. Martinez, P. Norberg, A. Fernandez-Soto, B. Ascaso, A. I. Merson, J. A. L. Aguerri, F. J. Gastander, 
L. Hurtado-Gil, G. Lopez-Sanjuan, A. Molino, A. Montero-Dorta, M. Stefanon, E. Alfaro, T. Aparicio-Villegas, N. Benitez, 
T. Broadhurst, J. Cabrera-Gaiio, J. Gepa, M. Cervino, D. Gristobal-Hornillos, A. del Olmo, R. M. Gonzalez-Delgado, 

G. Husillos, L. Infante, I. Marquez, J. Masegosa, M. Moles, J. Perea, M. Povic, F. Prada, and J. M. Quintana. The 
ALHAMBRA survey: evolution of galaxy clustering since $z \sim 1$. ArXiv e-prints, November 2013. 

[3] S. Basilakos, M. Plionis, K. Kovac, and N. Voglis. Large-scale structure in the HI Parkes All-Sky Survey: filling the voids 
with HI galaxies? MNRAS, 378:301-308, June 2007. 

[4] R. G. Bower, A. J. Benson, R. Malbon, J. G. Helly, C. S. Frenk, G. M. Baugh, S. Cole, and C. G. Lacey. Breaking the 
hierarchy of galaxy formation. MNRAS, 370:645-655, August 2006. 

[5] E. Branchini. Probing the Mass Distribution with IRAS Galaxies. ArXiv Astrophysics e-prints, October 2001. 

[6] A. L. Goil, J. A. Newman, M. G. Cooper, M. Davis, S. M. Faber, D. C. Koo, and C. N. A. Willmer. The DEEP2 Galaxy 
Redshift Survey: Glustering of Galaxies as a Function of Luminosity at z = 1. Astrophys. J. , 644:671-677, June 2006. 

[7] S. Cole, W. J. Percival, J. A. Peacock, P. Norberg, C. M. Baugh, C. S. Frenk, I. Baldry, J. Bland-Hawthorn, T. Bridges, 
R. Cannon, M. Colless, and C. Collins. The 2dF Galaxy Redshift Survey: power-spectrum analysis of the final data set 
and cosmological implications. MNRAS, 362:505-534, September 2005. 

[8] J. Gomparat, E. Julio, J.-P. Kneib, G. Schimd, H. Shan, T. Erben, O. Ilbert, J. Brownstein, A. Ealet, S. EscofRer, 
B. Moraes, N. Mostek, J. A. Newman, M. E. S. Pereira, F. Prada, D. J. Schlegel, D. P. Schneider, and C. H. Brandt. 
Stochastic bias of colour-selected BAG tracers by joint clustering-weak lensing analysis. MNRAS, 433:1146-1160, August 
2013. 

[9] S. Contreras, C. Baugh, P. Norberg, and N. Padilla. How robust are predictions of galaxy clustering? ArXiv e-prints, 
January 2013. 



15 


[10] N. Dalai, O. Dore, D. Huterer, and A. Shirokov. Imprints of primordial non-Gaussianities on large-scale structure: Scale- 
dependent bias and abundance of virialized objects. Phys. Rev. D , 77(12):123514, June 2008. 

[11] C. di Porto, L. Amendola, and E. Branchini. Growth factor and galaxy bias from future redshift surveys: a study on 
parametrizations. MNRAS, 419:985-997, January 2012. 

[12] G. Di Porto, L. Amendola, and E. Branchini. Simultaneous constraints on bias, normalization and growth index through 
power spectrum measurements. MNRAS, 423:L97-L101, June 2012. 

[13] G. Di Porto, E. Branchini, J. Bel, F. Marulli, M. Bolzonella, O. Cucciati, S. de la Torre, B. R. Granett, L. Guzzo, 
G. Marinoni, L. Moscardini, U. Abbas, C. Adami, S. Arnouts, D. Bottini, A. Gappi, J. Goupon, I. Davidzon, G. De Lucia, 
A. Fritz, P. Franzetti, M. Fumana, B. Garilli, O. Ilbert, A. lovino, J. Krywult, V. Le Brun, O. Le Fevre, D. Maccagni, 
K. Malek, H. J. McGracken, L. Paioro, M. Polletta, A. Polio, M. Scodeggio, L. A. M. Tasca, R. Tojeiro, D. Vergani, 
A. Zanichelli, A. Burden, A. Marchetti, D. Martizzi, Y. Mellier, R. G. Nichol, J. A. Peacock, W. J. Percival, M. Viel, 
M. Wolk, and G. Zamorani. The VIMOS Public Extragalactic Redshift Survey (VIPERS). Measuring nonlinear galaxy 
bias at z'0.8. ArXiv e-prints, June 2014. 

[14] D. J. Eisenstein, W. Hu, and M. Tegmark. Cosmic Complementarity: H 0 and Omega M from Combining Cosmic 
Microwave Background Experiments and Redshift Surveys. ApJL , 504:L57, September 1998. 

[15] D. J. Eisenstein, H.-J. Seo, and M. White. On the Robustness of the Acoustic Scale in the Low-Redshift Clustering of 
Matter. Astrophys. J. , 664:660-674, August 2007. 

[16] H. A. Feldman, N. Kaiser, and J. A. Peacock. Power-spectrum analysis of three-dimensional redshift surveys. Astrophys. 
J. , 426:23-37, May 1994. 

[17] J. N. Fry and E. Gaztanaga. Biasing and hierarchical statistics in large-scale structure. Ap. J., 413:447-452, August 1993. 

[18] E. Gaztanaga, P. Norberg, C. M. Baugh, and D. J. Croton. Statistical analysis of galaxy surveys - II. The three-point 
galaxy correlation function measured from the 2dFGRS. MNRAS, 364:620-634, December 2005. 

[19] T. Giannantonio, C. Porciani, J. Carron, A. Amara, and A. Pillepich. Constraining primordial non-Gaussianity with future 
galaxy surveys. MNRAS, 422:2854-2877, June 2012. 

[20] M. Grossi, L. Verde, C. Carbone, K. Dolag, E. Branchini, F. lannuzzi, S. Matarrese, and L. Moscardini. Large-scale 
non-Gaussian mass function and halo bias: tests on N-body simulations. MNRAS, 398:321-332, September 2009. 

[21] H. Hoekstra, L. van Waerbeke, M. D. Gladders, Y. Mellier, and H. K. C. Yee. Weak Lensing Study of Galaxy Biasing. 
Astrophys. J. , 577:604-614, October 2002. 

[22] E. Huff, A. E. Schulz, M. White, D. J. Schlegel, and M. S. Warren. Simulations of baryon oscillations. Astroparticle 
Physics, 26:351-366, January 2007. 

[23] E. Julio, J. Rhodes, A. Kiessling, J. E. Taylor, R. Massey, J. Berge, C. Schimd, J.-P. Kneib, and N. Scoville. COSMOS: 
Stochastic Bias from Measurements of Weak Lensing and Galaxy Clustering. Astrophys. J. , 750:37, May 2012. 

[24] 1. Kayo, Y. Suto, R. C. Nichol, J. Pan, 1. Szapudi, A. J. Connolly, J. Gardner, B. Jain, G. Kulkarni, T. Matsubara, R. Sheth, 
A. S. Szalay, and J. Brinkmann. Three-Point Correlation Functions of SDSS Galaxies in Redshift Space: Morphology, 
Color, and Luminosity Dependence. PASJ, 56:415-423, June 2004. 

[25] L. Knox, R. Scoccimarro, and S. Dodelson. Impact of Inhomogeneous Reionization on Cosmic Microwave Background 
Anisotropy. Physical Review Letters, 81:2004-2007, September 1998. 

[26] K. Kovac, C. Porciani, S. J. Lilly, C. Marinoni, L. Guzzo, O. Cucciati, G. Zamorani, A. lovino, P. Oesch, M. Bolzonella, 
Y. Peng, and B. Meneux. The Nonlinear Biasing of the zCOSMOS Galaxies up to z ' 1 from the 10k Sample. Astrophys. 
J. , 731:102, April 2011. 

[27] R. Laureijs, J. Amiaux, S. Ardulni, J. . Augueres, J. Brinchmann, R. Cole, M. Cropper, C. Dabin, L. Duvet, A. Ealet, and 
et al. Euclid Definition Study Report. ArXiv e-prints, October 2011. 

[28] A. Lewis, A. Challinor, and A. Lasenby. Efficient Computation of Cosmic Microwave Background Anisotropies in Closed 
Friedmann-Robertson-Walker Models. Astrophys. J. , 538:473-476, August 2000. 

[29] C. Marinoni, O. Le Fevre, B. Meneux, A. lovino, A. Polio, O. Ilbert, G. Zamorani, L. Guzzo, and A. Mazure. The VIMOS 
VLT Deep Survey. Evolution of the non-linear galaxy bias up to z = 1.5. A&A , 442:801-825, November 2005. 

[30] F. Marulli, M. Bolzonella, E. Branchini, 1. Davidzon, S. de la Torre, B. R. Granett, L. Guzzo, and A. lovino. The 
VIMOS Public Extragalactic Redshift Survey (VIPERS). Luminosity and stellar mass dependence of galaxy clustering at 
0.5<z<l.l. ArXiv e-prints, March 2013. 

[31] T. Nishimichi, 1. Kayo, C. Hikage, K. Yahata, A. Taruya, Y. P. Jing, R. K. Sheth, and Y. Suto. Bispectrum and 
Nonlinear Biasing of Galaxies: Perturbation Analysis, Numerical Simulation, and SDSS Galaxy Clustering. PASJ, 59:93- 
106, February 2007. 

[32] P. Norberg, C. M. Baugh, E. Hawkins, S. Maddox, D. Madgwick, O. Lahav, S. Cole, C. S. Frenk, 1. Baldry, J. Bland- 
Hawthorn, T. Bridges, R. Cannon, M. Colless, C. Collins, W. Couch, G. Dalton, R. De Propris, S. P. Driver, G. Efstathiou, 
R. S. Ellis, K. Glazebrook, C. Jackson, 1. Lewis, S. Lumsden, J. A. Peacock, B. A. Peterson, W. Sutherland, and K. Taylor. 
The 2dF Galaxy Redshift Survey: the dependence of galaxy clustering on luminosity and spectral type. MNRAS, 332:827- 
838, June 2002. 

[33] P. Norberg, C. M. Baugh, E. Hawkins, S. Maddox, J. A. Peacock, S. Cole, C. S. Frenk, and J. Bland-Hawthorn. The 2dF 
Galaxy Redshift Survey: luminosity dependence of galaxy clustering. MNRAS, 328:64-70, November 2001. 

[34] S. E. Nuza, A. G. Sanchez, F. Prada, A. Klypin, D. J. Schlegel, S. Gottloeber, A. D. Montero-Dorta, and M. Manera. The 
clustering of galaxies at z'0.5 in the SDSS-IH Data Release 9 BOSS-CMASS sample: a test for the LCDM cosmology. 
ArXiv e-prints, February 2012. 

[35] D. Schlegel, F. Abdalla, T. Abraham, C. Ahn, C. Allende Prieto, J. Annis, E. Aubourg, and M. Azzaro. The BigBOSS 
Experiment. ArXiv e-prints, June 2011. 



16 


[36] A. E. Schulz and M. White. Scale-dependent bias and the halo model. Astroparticle Physics, 25:172-177, March 2006. 

[37] U. Seljak. Redshift-space bias and 13 from the halo model. MNRAS, 325:1359-1364, August 2001. 

[38] H.-J. Seo and D. J. Eisenstein. Probing Dark Energy with Baryonic Acoustic Oscillations from Future Large Galaxy 
Redshift Surveys. Astrophys. J. , 598:720-740, December 2003. 

[39] H.-J. Seo and D. J. Eisenstein. Baryonic Acoustic Oscillations in Simulated Galaxy Redshift Surveys. Astrophys. J. , 
633:575-588, November 2005. 

[40] H.-J. Seo and D. J. Eisenstein. Improved Forecasts for the Baryon Acoustic Oscillations and Cosmological Distance Scale. 
Astrophys. J. , 665:14-24, August 2007. 

[41] P. Simon, M. Hetterscheidt, M. Schirmer, T. Erben, P. Schneider, C. Wolf, and K. Meisenheimer. GaBoDS: The Garching- 
Bonn Deep Survey. VI. Probing galaxy bias using weak gravitational lensing. A&A , 461:861-879, January 2007. 

[42] R. A. Skibba, M. S. M. Smith, A. L. Coil, J. Moustakas, J. Aird, M. R. Blanton, A. D. Bray, R. J. Cool, D. J. Eisenstein, 
A. J. Mendez, K. C. Wong, and G. Zhu. PRIMUS: Galaxy Clustering as a Function of Luminosity and Color at 0.2<z<l. 
ArXiv e-prints, October 2013. 

[43] R. E. Smith, R. Scoccimarro, and R. K. Sheth. Scale dependence of halo and galaxy bias: Effects in real space. Phys. Rev. 
D , 75(6):063512, March 2007. 

[44] V. Springel, S. D. M. White, A. Jenkins, C. S. Frenk, N. Yoshida, L. Gao, J. Navarro, R. Thacker, D. Croton, J. Helly, 
J. A. Peacock, S. Cole, P. Thomas, H. Couchman, A. Evrard, J. Colberg, and F. Pearce. Simulations of the formation, 
evolution and clustering of galaxies and quasars. Nature (London) , 435:629-636, June 2005. 

[45] M. E. C. Swanson, M. Tegmark, M. Blanton, and I. Zehavi. SDSS galaxy clustering: luminosity and colour dependence 
and stochasticity. MNRAS, 385:1635-1655, April 2008. 

[46] M. Tegmark. Measuring Cosmological Parameters with Galaxy Surveys. Physical Review Letters, 79:3806-3809, November 
1997. 

[47] M. Tegmark and B. C. Bromley. Observational Evidence for Stochastic Biasing. ApJL , 518:L69-L72, June 1999. 

[48] L. Verde, A. F. Heavens, W. J. Percival, S. Matarrese, C. M. Baugh, J. Bland-Hawthorn, and T. Bridges. The 2dF Galaxy 
Redshift Survey: the bias of galaxies and the density of the Universe. MNRAS, 335:432-440, September 2002. 

[49] I. Zehavi, Z. Zheng, D. H. Weinberg, J. A. Frieman, A. A. Berlind, M. R. Blanton, R. Scoccimarro, R. K. Sheth, M. A. 
Strauss, I. Kayo, Y. Suto, M. Fukugita, O. Nakamura, N. A. Bahcall, J. Brinkmann, J. E. Gunn, G. S. Hennessy, Z. Ivezic, 
G. R. Knapp, J. Loveday, A. Meiksin, D. J. Schlegel, D. P. Schneider, I. Szapudi, M. Tegmark, M. S. Vogeley, D. G. York, 
and SDSS Collaboration. The Luminosity and Color Dependence of the Galaxy Gorrelation Function. Astrophys. J. , 
630:1-27, September 2005. 

[50] The catalogs are publicly available at http://astro.dur.ac.uk/~d40qra/lightcones/EUGLID/ 


